Quantum escape of sudden future singularity 
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Abstract 

We consider finite-time, future (sudden or Big Rip type) singu- 
larities which may occur even when strong energy condition is not 
violated but equation of state parameter is time-dependent. Recently, 
example of such singularity has been presented by Barrow, we found 
another example of it. Taking into account back reaction of conformal 
quantum fields near singularity, it is shown explicitly that quantum 
effects may delay (or make milder) the singularity. It is argued that 
if the evolution to singularity is realistic, due to quantum effects the 
universe may end up in deSitter phase before scale factor blows up. 
This picture is generalized for braneworld where sudden singularity 
may occur on the brane with qualitatively similar conclusions. 
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1. Classical sudden singularities. The knowledge of the dark energy 
equation of state parameter w is extremely important due to various rea- 
sons. According to current astrophysical data this parameter is close (above 
or below) to —1. If it less than —1 then such dark energy (of phantom sort) 
drives the universe to finite time future singularity (Big Rip) [1, 2, 3] with 
catastrophic consequences for future civilizations. According to various esti- 
mations in this case the universe will exist for several more billions of years 
at best, before ending up in a cosmic doomsday. How realistic is such sudden 
future singularity? Is it common property of only dark energy with w less 
than —1? The partial answer to second question is given in the recent paper 
[4] where it was shown that even for usual matter (and when strong energy 
condition is not violated!) the sudden future singularity may occur. 

In the present Letter we argue that near to future singularity the quan- 
tum effects may become dominant and they can stop (or delay) the finite 
time future singularity (for general discussion of classical singularities, see 
book[5]). The same occurs for brane sudden singularities. This indicates 
that future Quantum Gravity era is quite possible and deserves some atten- 
tion. In principle, the proof that sudden singularity does not occur maybe 
quite simple. Indeed, one should write quantum-corrected FRW equations 
at late time and show that there is stable non-singular solution (at some 
conditions). For instance, such strategy in case of trace anomaly driven in- 
flation model shows that universe evolves to non-singular deSitter space. We 
go even further, taking into account quantum effects contribution as back- 
reaction near the future singularity. Even in such complicated formulation, 
there are indications that finite time singularity is soften (or delayed). 

Let us consider the generalization of Barrow's model [4]. In this model, 
the matter has been given implicitly via the FRW equations: 




(1) 



The spatially-flat FRW metric is considered 




(2) 



i=l,2,3 



and H = Eq.(l) shows that p is always positive. 
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We now assume H has the following form: 



H(t) = H(t) + A'\t s -t\ a . (3) 

Here H is a smooth, differentiable (infinite number of times, in principle) 
function and A' and t s are constants. Another assumption is that a constant 
a is not a positive integer. Then H(t) has a singularity at t — t s . In case 
a is negative integer, the singularity is pole. Even if a is positive, in case 
a is not an integer, there appears cut-like singularity. It is important for 
us that singularity presents. If we can consider the region t > t s , there is 
no finite-time future singularity. One also assumes the singularity belongs to 
the type (3) although there are other types of singularities, say, logarithmic 
ones. 

When a > 1, one gets 

p r^-p^ -H(t s f . (4) 

Hence w = - = — 1, which may correspond to the positive cosmo logical 
constant. 

The case < a < 1 corresponds to Barrow's model and when t ~ t s , we 
find 



6 ir i+\ 2 AA'a Q _! 

p~—H(t s ) , p~± — —\t s -t\ . (5) 



Here the plus sign in ± corresponds to t < t s case and the minus one to 
t > t s . In the following, the upper (lower) sign always corresponds to t < t s 
(t > t s ). The parameter of equation of state w is given by 

w = ± U J . 6 

3 H(t s f V ' 

Then it follows that w is positive in two cases: one is A' > and t < t s , 
which directly corresponds to Barrow's model, and another is A' < t > t s . 
In other cases, w is negative. 

When — 1 < a < 0, the energy density p and the pressure is given by 

P=—2-\ts-t\ , —\t s -t\ . (7) 
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The parameter of equation of state is 

w = ±^\t.-tr~ 1 , (8) 

which diverges at t — t s . Here w is positive when A' > and t < t s or 
A' < and t > t s . The former case corresponds to future singularity even 
if w is positive. The singularity can be regarded as a Big Rip. (For the 
recent comparison of phantom Big Rip with above type of it, see [6]. The 
mechanisms to escape of phantom Big Rip were suggested in refs.[7, 10].) 
The case a = — 1 gives 



6A 



,2 



K 2 



t s - if 2 , p ~ ~ (±2A' + \t s - t\~ 2 , (9) 



which may correspond to the scalar field with exponential potential. The 
parameter w is given by 

w = -l T ^. (10) 

Near t — t s , the universe is expanding if A' > and t < t s or A' < and 
t > t s . The former case corresponds to the phantom with w < —1. In the 
latter case, if 2 > A' > 0, the equation of state describes the usual matter 
with positive w and if A' > 2, the matter may be the quintessence with 
> w > -1. 

If a < —1, one obtains 



6A 



l2 



P= ~P= \t s ~ A 2 " , (11) 



K 2 



which gives w = —1 as for the cosmological constant case. In this case, 
however, there is rather strong singularity at t — t s . 

We should note that w is not constant in general. When w is constant, 
except the case w = — 1, the energy density p and the pressure p are given by 
(9) and always behave as \t s — 1|~ 2 . If we consider only expanding universe 
(not shrinking one), the future singularity appears when w < — 1, which 
corresponds to phantom. In the case with constant w > — 1, there does not 
appear the future singularity in the expanding universe although there is past 
or initial singularity at t — t s . 

Notice that usually w — — 1 corresponds to cosmological constant pres- 
ence but in general, in case that w — > — 1 when t — > t s there can appear a 
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singularity in the expanding universe as in (4) and (11). In case of (4), the 

f (or a, f , and £?; 



singularity is rather mild, and H and ^M- (or a, and Cf) are finite and 



(or ) diverges at t — t s . On the other hand, the case (11) shows very 

strong singularity, where the scale factor a behaves as a oc e + «+ 1 ' s 1 

On the contrary, when w is not constant but varies with t, even if the 
strong energy condition 

Prn > , p m + 3p m > (12) 

is satisfied, there can be a future singularity. An example is given in [4], 
which corresponds to (5). We have found another example given in (7) with 
A' > and t < t s . 

2. Quantum effects near Big Rip. When approaching the sudden fu- 
ture singularity, the curvature invariants grow and typical energies increase. 
This indicates that quantum effects (or even quantum gravity effects) should 
become significant near future singularity. They are usually associated with 
higher derivative terms. Imagine that quantum effects of conformal fields 
are dominant, as it occurs in some models of early universe, like in the 
trace anomaly driven inflation [8]. Let us take into account the conformal 
anomaly contribution as back-reaction near singularity. The correspondent 
energy density pa and pressure pa are (see, for instance, [9]) 

^-^-Q,-"){-^-^f + 3(f)1 (13) 




-12H- 2 " •"" ldH - (dH ^ 



dt 2 



-(f) H 



In general, with N scalar, Ni/ 2 spinor, Aq vector fields, N 2 (= or 1) gravi- 
tons and A^hd higher derivative conformal scalars, b, b' and b" are given by 

_N + 6N 1/2 + 12Aq + 611A^ 2 - 8A^ HD 
b ~ 120(4tt) 2 

Ar + llA^ 1/2 + 62Ar 1 + 14117V 2 -28A^HD u „ n , 1K , 
b= 36W ' (W) 
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We should note b > and b' < for the usual matter except the higher 
derivative conformal scalar s. 

Near the (future) singularity, the curvatures and their derivatives with 
respect to t become large. Since the quantum correction includes the powers 
of the curvatures and their t-derivatives, the correction becomes large and 
important near the singularity. It may occur that it delays (or stop) the 
singularity in the analogy with trace anomaly driven inflation [8] . 

As the first example, Barrow's model [4] is discussed. Following [4], we 
first consider the case that the scale factor a(t) is given by 

a(t) = A + Bt q + C(t s -t) n . (16) 

Here A > 0, B > 0, q > 0, and t s > are constants and C = —At~ n . One 
also assumes t < t s and 2 > n > 1. There appears a singularity when t — > t s , 
where — > +oo. Since 

a dt z 



qBt^ 1 - Cn (t 8 - tf' 1 qBtj- 1 Cn (t s - t)"" 1 
~ A + Bt<i + C (t 8 - t) n ~ A + Bti A + Bf s 

when t ~ t s , the parameters in (3) correspond to 



Cn 

A' = — — - , a = n - 1 . (18) 

A + Bti' 



The classical FRW equations 



6 2 2 rf2 ° Pm + 3Pm / in x 



/.■■ « 2 a <lt' 2 6 

show that 



6g2jB2t 2 ? - 2 Cn(n-l)(t s -tr 2 

> , p m ~ — q > , (20) 



Km K 2 (A + StI) ' ^ m A + Bti 

near the singularity £ ~ t s . The energy density p m is finite but the pressure 
p m diverges. Singularity presents but the strong energy condition (12) is 
satisfied since p and p are positive. 

With the account of the quantum corrections (13) and (14), the classical 
FRW equations (19) are modified as 



6 2 2 d 2 a p m + p A + 3p m + 3p A 

—tl = p m + Pa , -ir- 



2 H =p m + p A , — n ^i = g • ( 21 ) 
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As clear from the expressions (13) and (14), pA and pa become singular at 
t = t s . It follows from (21) 



V3 J (A + Btif 

Pm ~ -Pa 

~ 2^ + y) C "'"- 1 ><"- 2 ><"; 3) "-- <r4 <0. (22) 

Thus, the energy density p m is positive but diverges at t — t s . Since the pres- 
sure p m is negative and more singular than p m , the strong energy condition 
(12) is violated and also p + p < 0. 

For the example (20), the singularity at t — t s might become milder due 
to the quantum correction. To see this, we assume 4 > n > 3, instead of 
2 > n > 1, in (16). For this choice a, ^, , and are finite when t — t s 
but ^ diverges there. Then the singularity is rather milder than the case 
2 > n > 2, where ^ diverges. Using (21) around t = t s , one gets 

p m = —H - p A 

K 

6q 2 B 2 t 2 s «- 2 6b'BYt 4 s i- 4 



k 2 (A + Bti) (A + Btif 

-sf-b + b"^ ( B 2 q 2 ( q -l)( q -:n^- s 2B*<f\<i - 1W- 



.3 J\ (A + Bti) 2 (A + Bti) 3 

3B 4 qHi g -* 



(A + Bti) 4 
Pm ~ ~Pa 

2_ b \ Cn(n-l)(n-2)(n-3) (t. - t) n ~ 4 
3 J A + Bf s 

In the above expression, p m is not positive in general. However, in a special 
case b = —b' > 0, q — 1, it reduces to 

(24, 



k 2 (A + Bt s ) ' 
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which is positive. p m in (23) is, however, negative since C = —At~ n < in 
[4]. In principle, there is no any strong reason to exclude the case that C is 
positive. In the case C > 0, if one replaces 

2 + Cn(n- l)(n - 2)(ra - 3) -> -Cn(n - 1) , 

ra -> ra + 2 , (25) 

in p m (23), we obtain p m (20). The arguments may be reversed. If we start 
with the matter (20) (say, with q — 1), due to the quantum correction, we 
obtain the following scale factor a (instead of (16)): 

a(t) =A + Bt + C'(t a - tf . (26) 

Here A > 0, B > 0, C > 0, and t s > are constants and 4 > n' > 3. This 
demonstrates that the singularity becomes much milder due to the quantum 
correction. 

In general, near the singularity pA becomes dominant. The natural as- 
sumption is then 

P ~ ~Pa ■ (27) 

For the case given in (7), we find 

H{t) = H'{t) + A' \t s - t\ 2a+l , (28) 

Here H'{t) is a smooth and differentiable function and A' satisfies the follow- 
ing equation: 

^ = -2 06 + b'^j H' (t a ) (2a + 1) aA' > . (29) 

For the matter corresponding to (7), H behaves as H ~ |t s — but with 

the account of the quantum correction, one obtains H ~ |t s — tp 1 " 1 . Thus, 
the singularity becomes milder under quite general conditions. One should 
not forget that account of quantum gravity corrections (see recent exam- 
ple in [10]) may improve the situation with the singularity occurence even 
more. Finally, the quantum effective action which is usually non-local and 
higher- derivative typically leads to violation of energy conditions on quantum 
level. Hence, in many cases it may support the singularities. Nevertheless, 
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it is interesting that like in anomaly driven inflation there occur the sta- 
ble regimes where quantum effects may drive (at least, partially) the future 
universe out of the singularity Moreover, like in anomaly-driven inflation 
(with correspondent boundary conditions) one may argue that even future 
deSitter universe is quite realistic possibility as the corresponding solution of 
quantum-corrected FRW equation exists. 

3. Quantum escape in the brane-world. It is interesting that the 
above consideration of sudden future singularity may be easily generalized 
for brane-world models (for recent review, see [11, 12, 9] and refs therein). 
Let the 3-brane is embedded into the 5d bulk space as in [13]. Let g^ be 
the metric tensor of the bulk space and be the unit vector normal to the 
3-brane. The metric q^ induced on the brane has the following form: 

%v = - Jfyra„ . (30) 

The initial action is 

S = J d 5 x^g~ [^R {b) - 2A + • • • j + S branc (g) . (31) 

Here, the 5d quantities are denoted by the suffix ^5) and 4d ones by *4). In 
(31), ••• expresses the matter contribution and Sbrane is the action on the 
brane, which will be specified later. The bulk Einstein equation is given by 

^ (*{2 " = T m (32) 

If one chooses the metric near the brane as: 

ds 2 = dx 2 + q^dafdx" , (33) 

the energy-momentum tensor T^ u has the following form: 

T m = T^ lk mattcr - kg, v + 5( X ) (-Xq,u + r^) . (34) 

Here T^ lk mattcr is the energy-momentum tensor of the bulk matter, A is 
the bulk cosmological constant, A is the tension of the brane, and r^ u ex- 
presses the contribution due to brane matter. Without the bulk matter 
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^j-buik matter _ following the procedure in [13], the bulk Einstein equation 



- /IV 

can be mapped into the equation on the brane: 



= "2 [ A + "V J ^ + V + ^ - . (35) 

Here 7i^ u is given by 

= ~\v v a + + ^q^T aP T af5 - ^q^r 2 . (36) 

On the other hand, E^ u is defined by the bulk Weyl tensor C)Jpa but as the 
bulk spacetime is anti-deSitter one, where the Weyl tensor vanishes, we put 
E^y = 0. Like in [14] the tension of the brane A is fine-tuned to be 

= A + M^! . (37) 

The metric on the brane is supposed to be the FRW metric (2). Let the 
energy density p m and the pressure p m of the matter on the brane are given 
as in (20): 

Pm Pm 2r I a. a\— e /oo\ 

Ttt = — , Tij = —CL Oij , p m ~ po , p m ~ p [t s - t) . (38) 

Here p > 0, po > 0, and 1 > e > are constants. Then the strong energy 
condition (12) is satisfied. Now the effective Einstein equation (35) has the 
following form: 



6 tj2 , K 1 2 ^ d 2 a p m + 3p m k\ , 2 \ ; on \ 

-H = p m + - Pm , — — = _ + - (3 PmPm + 2 Pm ) . (39) 

Here 



k 2 ^^ . (40) 



When t ~ t s , one gets 




Pm . PO f . . K 2 5 po , , 
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which shows qualitatively the same behavior with that given by (16) with 
the following identification 



6q 2 BH 2 s q - 2 Cn(ra-l) 

Po = kha+bi*) ' n = -^vwr^ €=2 - n - (42) 



Hence, even for the brane- world, the qualitative behavior of the singularity is 
not changed (for earlier discussion of finite time singularities in brane-world, 
see [15]). 

The contribution from the conformal anomaly as a quantum correction 
from the matter on the brane may be again included. This can be done by 
replacing p m and p m in (39) with p m + p\ and p m + pa, respectively. Near 
singularity t ~ t s , one gets 

Pm ~ -Pa , (43) 



that is, 



which gives 



p„(t s -*r~2|'|6+f,")^f , (44) 



"={h b ')w^0w^ +m - (45) 

Here H(t) is a smooth, differentiable (at least, three times) function. Hence, 
the qualitative behavior is similar to that given by (26) with identification 

e = 4 - ri . (46) 

Thus, the singularity becomes much milder. The quantum effects of confor- 
mal matter help to realize the scenario of brane inflation called Brane New 
World [16, 17, 18]. As it follows from above consideration, the account of 
such quantum effects may help to avoid (or, at least to delay) the future 
sudden singularity. Moreover, as anomaly driven brane inflationary solution 
of FRW equations exists [16, 17, 18], one can argue that at similar conditions 
the future brane universe ends up in deSitter phase. 

The final picture (if realistic) looks quite strange. Quantum gravity era 
of the early universe helps in the realization of inflation. Subsequently, af- 
ter quite complicated history the universe evolves to dark energy dominated 
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phase which eventually may evolve to sudden singularity. The approach to 
finite-time future singularity brings back to life the Quantum Gravity era. In 
its turn, quantum effects may stop further evolution into the singularity what 
decreases the typical energies. In resulting, non-singular universe the quan- 
tum effects become negligible. Eventually, new circle of inflation (deflation) 
starts. This all reminds very much oscillating universe. 
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